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On the Diffraction of Light by Spheres of Small Relative 

Index. 

By Lord Eayleigh, O.M., F.E.S. 
(Eeceived February 4, — Eead February 26, 1914.) 

In a short paper " On the Diffraction of Light by Particles Comparable 
with the Wave-length/'* Keen and Porter describe curious observations upon 
the intensity and colour of the light transmitted through small particles of 
precipitated sulphur, while still in a state of suspension, when the size of the 
particles is comparable with, or decidedly larger than, the wave-length of 
the light. The particles principally concerned in their experiments appear 
to have decidedly exceeded those dealt with in a recent paper,f where the 
calculations were pushed only to the point where the circumference of the 
sphere is 2*25 X. The authors cited give as the size of the particles, when 
the intensity of the light passing through was a minimum, 6 to 10 //,, that is 
over 10 wave-lengths of yellow light, and they point out the desirability "of 
extending the theory to larger spheres. 

The calculations referred to related to the particular case where the 
(relative) refractive index of the spherical obstacles is 1*5. This value was 
chosen in order to bring out the peculiar polarisation phenomena observed in 
the diffracted light at angles in the neighbourhood of 90°, and as not inappro- 
priate to experiments upon particles of high index suspended in water. 
I remarked that the extension of the calculations to greater particles would 
be of interest, but that the arithmetical work would rapidly become heavy. 

There is, however, another particular case of a more tractable character^ 
viz., when the relative refractive index is small ; and although it may not be 

* ' Roy. Soc. Proc.,' A, vol. 89, p. 370 (1913). 

f < Roy. Soc. Proc.,' A, vol. 84, p. 25 (1910) ; < Scientific Papers,' vol. 5, p. 547. 
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the one we should prefer, its discussion is of interest and would be expected 
to throw some light upon the general course of the phenomenon. It has 
already been treated up to a certain point, both in the paper cited and the 
earlier one* in which experiments upon precipitated sulphur were first 
described. It is now proposed to develop the matter further. 

The specific inductive capacity of the general medium being unity, that of 
the sphere of radius E is supposed to be K, where K— 1 is very small. 
Denoting electric displacements by/, g, h, the primary wave is taken to be 

h = e int e ikx , (1) 

so that the direction of propagation is along x (negatively), and that of 
vibration parallel to z. The electric displacements (fi, gi, hi) in the scattered 
wave, so far as they depend upon the first power of (K — l), have at a great 
distance the values 

in which 

f (* r 

P= — (K — l).e int e} k ^ x ~ r )dxdydz. (3) 

In these equations r denotes the distance between the point (a, /3, 7) where 
the disturbance is required to be estimated, and the element of volume 
(dxdydz) of the obstacle. The centre of the sphere E will be taken as the 
origin of co-ordinates. It is evident that, so far as the secondary ray is 
concerned, P depends only upon the angle (y) which this ray makes with the 
primary ray. We will suppose that % = in the direction backwards along 
the primary ray, and that % = 7r along the primary ray continued. The 
integral in (3) may then be found in the form 

jfeoosig J Ji ( 2m cos 2 x • cos <i>) cos2 <i> d <$>> ( 4 ) 



p = -: 



r now denoting the distance of the point of observation from the centre of the 
sphere. Expanding the Bessel's function, we get 

4ttE 3 (K — 1) e*< tt *~* r > r 1 __ wL 4- m * m 6 

3 L 275 2.4.5.7 2.4.6.5.7.9 

77b 1 

+ 2.4. 6.8.5.7. 9. Ii~~'"J J (5) 

in which m is written for 2kU cos ^x- ^ ^ s to. be observed that in this 
solution there is no limitation upon the value of E if (K—l) 2 is neglected 
absolutely. In practice it will suffice that (K — 1)E/X be small, \ (equal to 
2 Trfk) being the wave-length. 

* 'Phil. Mag., 5 vol. 12, p. 81 (1881) ; * Scientific Papers, 5 vol. 1, p. 518. 
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These are the formulae previously given. I had not then noticed that the 
integral in (4) can be expressed in terms of circular functions. By a general 
theorem due to Hobson* 






K , ,s o , , , // 7r \ T . x smm cosm /nx 

Jx (m cos <£) cos 2 <^</> = V Urn/ t( m ) = -^ ™~' ( ^ 

so that 

p = _rK~l).47rE^^^~^)( /s ^~^!^i N ) j (7) 

\ m 3 m 2 / 

in agreement with (5). The secondary disturbance vanishes with P, viz., 
when tan m = m, or 

m = 2*Eoos J x = w (1-4303, 2'4590, 3'4709, 4-4747, 5*4818, etc.).f (8) 

The smallest value of &E for which P vanishes occurs when % = 0, i.e. in the 
direction backwards along the primary ray. In terms of X the diameter is 

2E = 0-715X (9) 

In directions nearly along the primary my forwards, cosJx is small, and 
evanescence of P requires much larger ratios of E to X. As was formerly 
fully discussed, the secondary disturbance vanishes, independently of P, in 
the direction of primary vibration (a = 0, /3 = 0). 

In general, the intensity of the secondary disturbance is given by 

yS.+rf+W-dS^lJj). (10) 

in which P denotes P with the factor e*<»*-* r ) omitted, and is a function of ^, 
the angle between the secondary ray and the axis of x. If we take polar 
co-ordinates (^, cf>) round the axis of x 9 

1 ""^T == l"-siii a x co8 2 0; (11) 

and the intensity at distance r and direction (^, </>) may be expressed in terms 
of these quantities. In order to find the effect upon the transmitted light, 
we have to integrate (10) over the whole surface of the sphere r. Thus 



/yiJ 



rr f 71 " /& 2 p \ 2 

JJ sin X^X#(/i 2 +#i 2 + Ai 2 ) = Trjsinx^— °j (l + cos 2 % ) 

74/tr i\2ufif 7r - 7 /i , 2 N (sin m — m cos m) 2 
= irifc* (K - 1 ) 2 E 6 sm x d X (1 + cos 2 *) v ~ Q } ~ 

Jo ' fb 

_ 2 tt/,E (K-l) J^ w6 |2 j£ _ + __j. 

x {l+m 2 -f (m 2 — 1) cos 2m— 2m sin 2m}. (12) 

* « Lond. Math. Soc. Proc., 5 vol. 25, p. 71 (1893). 
t See < Theory of Sound,' § 207. 
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The integral may be expressed by means of functions regarded as known 
Thus on integration by parts 

{l + m 2 + (m 2 ~l) cos 2m — 2m sin 2m} 



o • 



m 



1 — cos 2 m sin 2 m __ 1 1 
4m 4 2m 3 2m 2 2' 



*m 







il-f-m 2 + (wi 2 — 1) cos 2m— -2m sin 2m} — 5 



+ 



2 m 2 . 



' ro l- cos 2m 7 cos 2m sin 2m ., 

AWM - r . 1 -1, 



m ' 2 m 2 m 



{l + m 2 + (m 2 — l)cos 2m — 2m sin 2m} 







m 



"ni 



1— cos 2m 7 t m 2 m sin 2m 5 cos 2m 5 



Accordingly, if m now stand for 2&B, we get 
sin % % ^<£ (/1 2 +^i 2 + hi 2 ) 



r 



2 ' 



= ^R 2 (K-1) 2 f _ 7 (1 -cos 2m) _ sin 2m + g + m2 

8 \ 2m 2 m 

If m is small, the { } in (13) reduces to 

+ Oxm 2 +^-m 4 , 

so that ultimately 

(13) = ^-tt/^E^K-I) 2 , (14) 

in agreement with the result which may be obtained more simply from (5). 
If we include another term, we get 

(13) = ^ 7r ^E6(K-l) 2 (l-^-\ (15) 

As regards the definite integral, still written as such, in (13), we have 
^Thzm^d^n = p|jL_g + J-...J dx= 7 + log(2m)-Ci(2m), 

° m ° (16) 

where 7 is Euler's constant (0*5772156) and Ci is the cosine-integral, 
defined by 

Ci(x)=p^d u . (17) 
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As in (16), when x is moderate, we may use 



xr 



Ci(a?) = 7 4- log (T— \ - — ~ -\-\ 



x q 



1.2 * 1.2.3.4 



(18) 



which is always convergent. When x is great, we have the semi-convergent 

series 

1 1.2 , 1.2.3.4 



Ci (x) = sin a? < ~ + 






x- 



X. X. . JL % O . X.Z.O.TC.O 

— COS& < ~ : H 






OF 



(19) 



Fairly complete tables of Ci (x), as well as of related integrals, have been 
given by Glaisher.* 

When m is large, Ci(2m) tends to vanish, so that ultimately 

* m l — cos2m 7 , , n x 

— dm — <y + log(2m). 



o 



m 



Hence, when kB, is large, (13) tends to the form 

(13) = |7rFE 4 (K-l) 2 . (20) 

Glaisher's Table XII gives the maxima and minima values of the cosine- 
integral, which occur when the argument is an odd multiple of |7r. Thus — 



n. 


Ci(mr/2). 


71. 


Ci(wr/2). 


1 
3 

5 


+ 0-4720007 
-0-1984076 
+ 0-1237723 


7 

9 

11 


-0-0895640 
+ 0-0700653 
-0 -05750 LI 



These values allow us to calculate the { } in (13), viz., 

4 



7(1 — cos 2m) sin 2m , K , « , /4 ,\ r f1 n ^. /n X1 /M . 
2^ '--^— +5+m 2 -f ^_~4j[ 7 + log 2m~Ci(2m)], (21) 

when 2m = n7r/2, and n is an odd integer. In this case cos 2m = and 
.sin 2m = ±1, so that (21) reduces to 



56 4 n 2 7T 2 I P>4- 

± ±+ l 5 + ^ + f4V-4 [7 + log(w7r/2)-.Ci(W2)]. (22) 



?i 2 7r 2 mr ' 16 * W 2 7r 2 



We find 



n. 


(22). 


n. 


(22). 


1 
3 
5 


-0530 
2-718 
10 -534 


7 

9 

11 


23 -440 
42 '382 
65 -958 



* ' Phil. Trans.,' vol. 160, p. 367 (1870). 
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For values of n much greater, (22) is sufficiently represented by ^ 2 7r 2 /16, 
or m 2 , simply. It appears that there is no tendency to a falling- off in the 
scattering, such as would allow an increased transmission. 

In order to make sure that the special choice of values for m has not 
masked a periodicity, I have calculated also the results when n is even. 
Here sin 2m = and cos 2m = + 1, so that (21) reduces to 



112 (lorO) iKi ^ 2 7T 2 , / 64 
-5—^ 4-0-1--— - + 



n 2 7r 2 



n 2 7r 2 



•4) [7 + log (*wr/2)-Ci (twt/2)]. (23) 



The following are required :- 



n. 


Ci (»»/2). 


». 


Ci(»ir/2). 


2 
4 
6 


+ -0738 
-0*0224 
+ -0106 


8 
10 


-0-0061 
+ 0*0040 



of which the first is obtained by interpolation from G-laisher's Table VI, an 
the remainder directly from (19). Thus : — 



n. 


(23). 


n. 


(23). 


2 

4 
6 


-7097 
6 -1077 
16 -156 


8 
10 


32 *336 
53 -477 



The better to exhibit the course of the calculation, the actual values of the 
several terms of (23) when n = 10 may be given. We have 

112 



u 2 tt 2 



= -0*11348, 
64 



7b 2 IT 2 

"16" 



61-685, 



^ 2 7T 2 



4-0*06485 = 3*93515, 



7 + l og (7r/2) + logn-Ci(^7r/2) = 0*57722 + 0*45158 + 2*30259 -0-0040 

= 13*094, 



so that 



64 



n 2 7r 2 



{7 + log (wtt / '2) -Ci(mr/ '2)} = 13*094. 



It will be seen that from this onwards the term n 2 7r 2 /l6, viz., m 2 , greatly 
preponderates ; and this is the term which leads to the limiting form (20). 

The values of 2E/X concerned in the above are very moderate. Thus> 
n = 10, making m = 47rE/X = IOtt/4, gives 2K/A, = 5/4 only. Neither 
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below this point, nor beyond it, is there anything but a steady rise in the 
value of (13) as X diminishes when E is constant. A fortiori is this the 
case when E increases and \ is constant. 

An increase in the light scattered from a single spherical particle implies, 
of course, a decrease in the light directly transmitted through a suspension 
containing a given number of particles in the cubic centimetre. The 
calculation is detailed in my paper " On the Transmission of Light through 
an Atmosphere containing Small Particles in Suspension,"* and need not be 
repeated. It will be seen that no explanation is here arrived at of the 
augmentation of transparency at a certain stage observed by Keen and 
Porter. The discrepancy may perhaps be attributed to the fundamental 
supposition of the present paper, that the relative index is very small, 
a supposition not realised when sulphur and water are in question. But I 
confess that I should not have expected so wide a difference, and, indeed, 
the occurrence of anything special at so great diameters as 10 wave-lengths 
is surprising. 

One other matter may be alluded to. It is not clear from the description 
that the light observed was truly transmitted in the technical sense. This 
light was much attenuated — down to only 5 per cent. Is it certain that it 
contained no sensible component of scattered light, but slightly diverted 
from its original course ? If such admixture occurred, the question would 
be much complicated. 

* ' Phil. Mag.,' vol. 47, p. 375 (1899) ; < Scientific Papers/ vol. 4, p. 397. 



